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ABSTRACT 

Optically very faint {R > 25.5) sources detected by the Spitzer Space Telescope at 
24/im represent a very interesting population at redshift z ~ [1.5 — 3]. They exhibit 
strong clustering properties, implying that they are hosted by very massive halos, 
and their mid-IR emission could be powered by either dust-enshrouded star-formation 
and/or by an obscured AGN. We report observations carried out with the MAMBO 
array at the IRAM 30m antenna on Pico Veleta of a candidate protocluster with five 
optically obscured sources selected from the 24/im Spitzer sample of the First Look 
Survey. Interestingly, these sources appear to lie on a high density filament aligned with 
the two radio jets of an AGN. Four out of five of the observed sources were detected. 
We combine these measurements with optical, infrared and radio observations to probe 
the nature of the candidate protocluster members. Our preliminary conclusions can be 
summarized as follows: the Spectral Energy Distributions of all sources include both 
AGN and starburst contributions; the AGN contribution to the bolometric luminosities 
ranges between 14 and 26% of the total. Such a contribution is enough for the AGN to 
dominate the emission at 5.8, 8 and 24/im, while the stellar component, inferred from 
SED fitting, prevails at 1.25mm and at A < 4.5/xm. The present analysis suggests a 
coherent interplay at high-z between extended radio activity and the development of 
filamentary large-scale structures. 

Key words: Galaxies: evolution, high redshift, ISM, clusters: general - ISM: general 



1 INTRODUCTION 

Optically obscured 24/im Spitzer sources are among the 
most strongly clustered objects in the Universe. Maglioc- 
chetti et al. (2007) studied a complete sample of 793 
optically faint {R > 25.5) sources with -F24nm ^ 0.35 mjy, 
drawn from the Spitz er Spac e Telescope First Look Sur- 
vey (FLS) (jPadda et al. 20061 '). Their two-point correlation 
function has an amplitude A = (7±2)-10"^ [w{e) = ylS""'*, 
with 9 in degrees], which is about a factor of 8 larger than 
that obtained for the w hole Spitzer F24,um ^ .35 mJy 
sample (see Fig. 6 of iMagliocchetti et al. 2007h . Their 
optical/mid-IR colours indicated that they are ultralumi- 
nous far- infrared galaxies at z ~ 1.6 — 2.7, a conclusion 
supported by redshift determinations based on infrared 
spectroscopy (IRS; e.g., Yan et al. 2005, 2007; Houck et 
al. 2005; Weedman et al. 2006a, Desai et al., 2007) and 



by comparisons with models for galaxy formation and 
evolution (e.g., iGranato et al. 20041 ). Their 3D comoving 
clustering radius was estimated to be ro = [15.2l2'g] Mpc 
[^(r) = (r/ro)"^'*] (Magliocchetti et al. 2007). 

Clustering results obtained on Bootes and UKIDSS 
fields by Brodwin et al. (2008) and Magliocchetti et al. 
(2008b) based on spectroscopic and photometric N(z) con- 
verge, within the errors, at giving similar high clustering 
lengths and halo masses (ro about 13 Mpc and masses 10^^'^- 
10^^ Mq for S>0.4 mjyjj. These measured clustering prop- 
erties indicate that these sources are associated to very mas- 



Note that the clustering length given by Magliocchetti et al. 
(2008b) is in Mpc for h = 0.7, not in Mpc, as stated by 

Brodwin et al. (2008) 
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sive dark matter (DM) halos, consistent with them inhab- 
iting the progenitors of groups/clusters of galaxies. Thus, 
these objects are important tools to investigate the abun- 
dance of very massive halos at high z and the evolution of 
large scale structure. 

Bright, high-z QSOs, powered by the most massive 
black-holes (BHs ), display similar clustering properties 
(IShen et al. 20071 ) and should therefore inhabit similar ha- 
los at z ~ 2. Obscured 24/im sources in the environment 
of bright QSOs (and/or radiogalaxies) at 2; ~ 2 then likely 
trace the largest density peaks and thus the very most mas- 
sive DM halos. This could be materialized by several fila- 
ments crossing at their position. High-z radiogalaxies proved 
to be particularly effective tracers of distant large scale over- 
densities that may collapse to form present-day massive clus- 
ters of galaxies (Pentericci et al. 1997; Venemans et al. 2007; 
Overzier et al. 2007; Miley & De Breuck 2008). Since both 
optically obscured 24/im Spitzer sources and high-z radio- 
galaxies are preferentially found within the rare large-scale, 
large-amplitude peaks in the primordial density fluctuation 
field, it is interesting to investigate possible links between 
these two source populations. 

A by-product of the clustering analysis by Magliocchetti 
et al. (2007) was the discovery of several galaxy alignments, 
reminiscent of the high-z filaments found in numerical simu- 
lations which describe the evolution of large-scale structure. 
In some cases, the candidate filaments were aligned with ra- 
dio jets, consistent with West's (1994) prediction that the 
radio axes reflect the developing large-scale clustering pat- 
tern in the early Universe. If so, these systems would be an 
important probe of the cosmological density fleld at early 
epochs. 

In this paper we start investigating the above issue by 
collecting all the available data for one of these candidate 
fllaments, and supplementing it with our own observations 
with the MAMBO array at the focus of the IRAM 30m 
antenna at 1.2mm. 

The addition of mm data is crucial in constraining the 
nature of the obscured 24/im sources. In fact their nature is 
currently under debate: a large fraction of them has radio 
spectral indices typi cal oi frustrated radio- loud A ctive Galac- 
tic Nuclei (AGNs) (|Maghocchetti et al. 2008^ 1 but their q 
values - i.e. the values of the ratio between far-IR (FIR) 
and radio fl ux densities - are similar to those of star forming 
objects (i.e. Ilbar et al. 20081 ). Their colours indicate mixed 
AGN and starburst properties (Magliocchetti et al. 2007; 
Pope et al 2008; Lonsdale et al. 2009). 

FIR/mm observations are very effective for determining 
the bolometric output by measuring the redshifted peak of 
the dust emission, for characterizing the broad band spec- 
trum which allows to separate AGN from starburst emis- 
sions, and for studying the parallel dust obscured evolution 
of super-massive BHs (SMBHs) and of their host galax- 
ies. The AGN contribution produces hotter dust and much 
higher F(mid-IR)/f (sub-mm) ratios than (sub)mm selected 
sources, which are generally high-z starburst galaxies. Lutz 
et al. (2005) found a mean ratio -F'i.2mm/i^24nm = 0.15 ± 
0.03 for a sample of AGN-dominated optically obscured 
Spitzer FLS sources with -F24/jm > 1 mjy. Spitzer obser- 
vations of starburst-dominated objects originally detected 
with SCUBA and MAMBO (Pope et al., 2006, Ivison et al.. 



2007, Lonsdale et al., 2008) and with F24^-m > 20/iJy, instead 
have a median Jsso^im /-F24pm 25 and -Fi.2mm/-F24Mm ^ 12. 

Estimating the AGN and star-forming components 
helps in investigating any evolutionary connection between 
them, as the relationship between the star-formation his- 
tory and the growth of the active nucleus is a key in- 
gredient for understanding the phys ics of the evolution 
of both massive galax ies and AGNs (|Granato et al. 20041 : 
ISchawinski et al. 20061 ). 

This work deals with a candidate protocluster system of 
optically obscured 24/im Spitzer FLS sources in the fleld of 
an extended radio-loud AGN. Such a system: i) provides 
us with a sample of sources with a 24^m flux distribu- 
tion which is representat ive of that found for the whole 
(|Magliocchetti et al. 20071 ) dataset; ii) allows us to inves- 
tigate in a direct way the transition between 'pure' star- 
forming galaxies and AGN-dominated sources as envisaged 
e.g. by the Granato et al. model; iii) perhaps most impor- 
tantly, it provides an excellent laboratory to study the large- 
scale interplay between AGN activity (such as that produced 
by radio jets) and star formation in the sources belonging 
to the protocluster. 

The central source of this complex radio sys- 
tem is identifled with FLSVLA172010. 2-^592425. 2 
(jGondon et al. 20031 ) and at 8" from it there is an- 
other radio source, 7C1719-I-5927 (|Hales et al 19951 ). The 
flux densities integrated over the whole system at 38, 330, 
610 MHz, and 1.4GHz are of 2000, 143, 33.7, and 35.6 
mJy, respectively. The 38, 330, and 1.4 GHz measurements 
are consistent with a power-law spectrum [S oc , with 
a ~ 1), while the 610 MHz flux density is a factor ~ 2 
below that. The central source has an optical counterpart 
with a magnitude oi R ^ 23.67 (Fadda et al. 2004). The 
spectroscopic redshift of this source is unknown. 

Section 2 reports all the data available on this sample 
and the new MAMBO observations. Results of our analysis 
are given in § 3 and discussed in § 4, while in § 5 we summa- 
rize our main conclusions. 



2 DATA 

The source sample used in this work was drawn from a par- 
ent sample whose selection criteria are outlined in this sec- 
tion. 



2.1 Sample selection and Spitzer data 

Of the original 793 R > 25.5 sources in the Maglioc- 
chetti et al. (2007) sample, 510 have IRAC photometry. 
On the basis of these data, sources were classifled ac- 
cording to their J8;jm/J24 Mm and F4.5^m/-F3.6Mm colours 
(jMagliocchetti et al. 2008ah as follows: 

- AGNO: F24Mm > 0.6 mJy; F8Mm/^24Mm > 0.1; 

- SBO (starburst): F24Mm ^ 0.6 mJy; F8^m/-F24pm < 0.1; 

- AGNl: F24Mm < 0.6 mJy; F^.si^m/ Fa.e^m > 4.5/3.6; 

- Unclassifled: ^24^111 < 0.6 mJy; F^.^t^m/ FsMt^m < 4.5/3.6 
or no IRAC information. 

The above classification is consistent with the results of mid- 
infrared (IRS) spectroscopy of samples of heavily obscured 
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Figure 1. Spitzer/YhA field centered at 17:20:05, +59:24:30. Filled circles correspond to the FLS Spitzer 24^m sources we observed 
with MAMBO. 



24^m-bright objects (brighter than 0.7-1 mjy). The ma- 
jority of them have spectral shapes indicative of a buried 
AGN (Martinez-Sansigre et al. 2006; Yan et al. 2005, 2007; 
Brand et al. 2007; Weedman et al. 2006a,b; Pope et al. 2008). 

2.2 Radio observations: selection of protoclusters 

Magliocchetti et al. (2008a) searched for radio counterparts 
to optically obscured Spitzer sources by cross-correlating the 
source sample of Magliocchetti et al. (2007) with the radio 
source catalogs by Condon et al. (2003) at 1.4GHz, Morganti 
et al. (2004) at 1.4GHz and Garn et al. (2007) at 610 MHz. 

A number of the selected sources turned out to lie in 
overdense regions of the sky in the neighborhood of radio- 
loud AGNs (see, e.g., Fig.[T}. To pursue our research further 
we then selected small complete subsamples in those high 
density regions with deep radio imaging and hosting a radio- 
loud AGN. In this paper we concentrate on the field shown in 
Fig. [T] The available Spitzer, radio and optical photometric 
data are listed in Table [1] 

2.3 Redshift constraints: IRAC and optical 
observations 

IRAC photometry can also be used to provide rough redshift 
estimates for those sources whose SEDs present a contribu- 



tion from an evolved stellar population in the 3.6/im and 
4.5/im bands. Based on th e form of the SEP in this wave- 
length range (see Fig. 6 in [Magliocchetti et al. 2008ah . and 
in particular on its interpolated slope (3 {F (x \^) between 
3.6/im and 4.5/im, we estimated that the redshifts of these 
objects are around z ~ 1.7. This makes it plausible that they 
are physically associated and not just close to one another 
due to projection effects. 

An independent assessment of the redshift range for this 
source population comes from the measured spectroscopic 
redshifts of optically faint 24/nm Spitzer sources by Dey et al. 
(2008) using the IRS Spitzer spectrograph and Keck optical 
spectroscopy. If we adopt Dey et al.'s colour criterion, our 
sources have a median redshift < z >— 1.8-1.9, roughly 
consistent with our previous estimate. 



2.4 MAMBO observations of the candidate 
protocluster 

Data were taken with the 117-channel Max Planck Millime- 
tre Bolometer Array (MAMBO; Kreysa et al. 1998) at the 
focus of the 30m IRAM antenna at Pico Veleta. Chopping 
photometry at 2 Hz with 32" throw in azimuth was gath- 
ered at 250 GHz (1.25 mm) in observing pool during Win- 
ter 2008 towards the sources listed in Tables [T]and (2] Data 
were obtained in good observing conditions (250 GHz atmo- 
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Source name 




3.6/xm 


4.5/im 


5.8/im 




24/jm 


1.4GHz 


610MHz 


ID 














^Jy 


piJy 


MJy 





J172037.75+592359.3 


>25.5 


29.93±4.33 


45.84±6.24 


<100 


<100 


400±80 


<100 


<182 


SB 


J172030.92+592308.3 


>25.5 


64.68±7.84 


53.51±7.15 


<100 


191.4±24.3 


390±75 


<100 


<182 


AGNO 


J172005.01+592430.0 


>25.5 


24.36±3.65 


33.52±5.18 


<100 


<100 


440±88 






AGNl 


J171939.17+592512.4 


>25.5 


20.89±3.37 


36.59±5.44 


<100 


<100 


380±74 


157±27 


<182 


AGNl 


J172025. 20+592308.0 


>25.5 


<100 


<100 


<100 


<100 


350±70 


<100 


<182 


UNCL 



Table 1. Optieally obseured 24/jm Spitzer sources observed with MAMBO during Winter 2008. The table lists the optical and Spitzer 
photometry available for these sources. Last two columns report the radio measurements when available. The source J172005+592430 
lies on top of a lobe of the radio source (see Fig. [TJ, that hides its possible radio emission. Upper limits are at the 3(7 level. The last 
column shows the type of identification assigned by Magliocchetti et al. (2007) to these sources. We have assigned in our analysis an 
error- bar of 20% to the 24/im fluxes (cfr text). 



Source name weighted first atmospheric average on source second atmospheric average on source 

flux observation transmission atm noise int time observation transmission atm noise int time 
(mjy) (mJy) t (mjys^^^) (min) (mJy) t (mjys^^^) (min) 



J172037.75+592359.3 


1.4±0.5 


1.2±0.5 


0.35 


36 


160 


3.4±1.4 


0.15 


90 


52 


J172030.92+592308.3 


<1.2 


0.12±0.63 


0.16 


50 


48 


-0.13±0.58 


0.36 


40 


64 


J172005.01+592430.0 


1.6±0.3 


1.0±0.5 


0.37 


46 


80 


2.2±0.5 


0.27 


63 


72 


J171939.17+592512.4 


1.4±0.5 


1.3±0.6 


0.5 


90 


102 


1.3±0.6 


0.3 


40 


52 


J172025.20+592308.0 


1.8±0.4t 


0.7±0.6 


0.35 


40 


60 


2.5±0.6 


0.26 


60 


80 



Table 2. MAMBO observations of 24 /im optically obscured Spitzer sources. The sources observed were previously classified - according 
to our colour selection criterion (Magliocchetti et al., 2007) — as either starbursts or obscured AGNs. The only undetected source was 
indeed classified as an AGN and its flux upper limit is 1.2 mJy (3cr). Atmospheric values in columns 4,5,8,9 refer to average values during 
the observations. Errors on the fluxes are statistical. Final flux is computed as a weighted average. 

f Because of the discrepancy between the two measurements of this source, we assign in our analysis an additional systematic error-bar 
of 20% to its mm flux. 



spheric opacity better than 0.2 and low sky noise), as re- 
sulting from the regular checks of atmospheric transmission 
and noise. Pointing was checked every hour with a nearby 
quasar. Each source has been observed in at least two inde- 

J iendent run s. Data were reduced with the MOPSIC software 
Zvlka 19981 ). 

To check for any malfunctioning of the detecting system 
with consequent instabilities of the bolometer output, we 
inspected subscan by subscan and reduced the data sepa- 
rately for each observing chunk. Subscans containing spikes 
(less than 10% of the total) were discarded. Observing logs 
with total integration times on sources, atmospheric param- 
eters, fluxes at each run and the final calibrated fluxes with 
statistical errorbars are reported in Tabled Longer integra- 
tion times were necessary because of the higher sky noise 
during some of the observations. Statistically, the sky noise 
was behaving properly as it scales with the square root of 
the integration time. 



3 RESULTS 

3.1 Interpreting the observed spectral energy 
distributions 

Figure [5] shows the photometric points from the R magni- 
tude down to the radio fluxes at 1.4 GHz and 610 MHz for 
the sources of our sample. This also includes the MAMBO 
observations presented in Section 2.4. The original 24^m er- 
rorbars given by Fadda et al (2006) for these sources are un- 
realistically small (1.5%); on the suggestion by the referee, 
we have adopted 20% error bars, as estimated for similar 
(80s) integration times. 

In Table |3] we report the results from two different types 
of fits and the corresponding values and luminosities. The 
first template used is the Arp220 SED, while the second has 
two components: Arp220 SED and an obscured torus tem- 
plate made available by PoUetta et al. (2007). In this latter 
case we could estimate the contribution to the bolometric lu- 
minosity of the two components. We find that the Arp 220 
Spectral Energy Distribution (SED) provides the best single 
component fit for all the observed SEDs. In fact, all of the 
30 SEDs of nearby star-forming galaxies presented by Vega 
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Figure 2. Spectral energy distributions of the four sources detected with MAMBO in our sample. Photometric points from the R 
magnitude down to the radio fluxes at 1.4 GHz and the 610 MHz are shown together with two spectral energy distributions which best 
fit the observed points: a starburst component represented by Arp220 SED (green dashed line) and an obscured torus (PoUetta ot al., 
2007, red dotted line). The SEDs arc shifted to a redshift of 1.73. The corresponding bolometric luminosities are 5-7 xlO^^ Lq (see 
Tabled. 



et al. (2008) and the other template SEDs made available 
by PoUetta (PoUetta et al. 2007) give substantially worse 
fits. By then using the Arp 220 SED as 'bona fide' tem- 
plate, we get a minimum fit at the same value of the 
redshift, z — 1.73, for each object belonging to the filament 
separately, with the exception of J172030. 9-1-592310 which 
is classified as an AGNO from IRAC photometry (see Table 
1). At this redshift, the sources have bolometric luminosities 
in the range 5-7-10^^ L0 and star formation rates of ~ 1000 
Mo/yr (see Table [sfl. Obviously, this surmise can only be 
confirmed by spectroscopic redshift measurements. 

No really good fit could be found for J172030.9-I-592310. 
In this case, neither the Arp 220 template, nor any other 
template mentioned above, can be accepted as a realistic de- 

^ The star formation rate was estimated from the total far-IR 
luminosity, which is essentially equal to the bolometric luminosity, 
using eq. (4) of Konnicutt (1998). 



scription of its SED. The lowest x^, still quite unsatisfactory 
(X^ = 23.9), fit was obtained by adding to the SED of Arp 
220 the PoUetta's "torus" template (heavily obscured AGN) 
at z = 1.214 (green line in Fig. [S]). Additionally, Figure [3] 
shows as a red line the best-fit with the same combination 
of SEDs red-shifted to z = 1.73. This shift in z substantially 
worsens the value (x^ ~ 48.5). 

In Table [3] we give the estimated bolometric luminosi- 
ties af the sources for the two models we have considered, i.e. 
pure starburst SED (I fit) and starburst plus AGN (II fit). 
The errors on the total luminosities were computed finding 
the contours in the parameter space for which x^ ~ Xmin + 1 
and computing the maximum and minimum luminosities on 
these contours. The resulting errors are of 30-35%. The main 
contribution to them comes from the uncertainties on red- 
shifts of the sources, which turn out to be Az ~0. 17-0.20. In 
the case of the SB-I-AGN fit we give the estimated bolomet- 
ric luminosities of each component. The global bolometric 
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Figure 3. Spectral energy distribution as in Figure |2] for the 
AGNO. The green dashed line shows the best fit, obtained adding 
the SED of Arp 220 with PoIIetta's "obscured torus" template, 
red-shifted to z = 1.214; the red solid line shows the best-fit 
(giving a much worse x^) derived with the same combination of 
SEDs red-shifted to z = 1.73. We classify this source as a type 2 
AGN candidate, i.e. an obscured AGN. 

luminosity is therefore the sum of the two contributions, and 
it is used to compute, in a consistent way, the SB and AGN 
fractions. The estimate of the bolometric luminosity clearly 
depends on the used SED template. No surprise that 2 sub- 
stantially different templates yield different results. Note, 
however, that differences are not very large. In the worst 
case (source J171939.1, leaving aside the pathological second 
source), going from a pure Arp220 template to a SB+AGN 
template decreases the bolometric luminosity by 36%. We 
do not attach any profound implication to the fact that for 
3 sources the bolometric luminosity decreases from case I to 
case II, while only in one case (last source) it increases: this 
is most likely happening by chance. 

Based on the SED analysis which is listed in Table 3, 
we find that the nature of our candidate protocluster mem- 
bers can be summarized as follows: the AGN contribution 
to the bolometric luminosities ranges between 14 and 26% 
of the total. A ~ 14% contribution to the total bolometric 
luminosity is enough for the AGN to dominate the emission 
at 5.8, 8 and 24/im, while the SB dominates at 1.25mm and 
at A < 4.5/im. The above fits are consistent with all these 
sources being at the same redshift. 

3.2 Nature of the obscured Spitzer 24fim sources 

The information provided by the spectral energy distribu- 
tion is not, however, compelling and the nature of these 
sources cannot be entirely evinced. However colour-colour 
plots may be further exploited to support our previous find- 
ings on the origin of the source emission. We compare in 
this section the colours of the sources presented in this work 
with those of other 24^m Spitzer sources, extracted with 




F(8/2m)/F(4.5jU,m) 



Figure 4. The colour— colour F{24fj,m) / F{8fj,m) vs. 
F{8fim) / F{A.5fim) plot is used as a diagnostic tool to discrimi- 
nate between starburst and AGN dominated SEDs. The vertical 
dashed line corresponds to the ratio (F(8/^m)/i<'(4.5/im) = 2, 
proposed by Pope et al. (2008) as the boundary between SB 
(having ratios < 2) and AGN dominated SEDs. Filled squares 
correspond to the Spitzer sources presented in this work (we 
note that for all but one case the 8fim measurement is a mere 
upper limit); filled triangles: Ivison et al. (2004); filled pentagon: 
Younger et al. (2007); filled circles: Sawicki & Webb (2005); 
crosses: Pope at. (2006); four-points stars: Lonsdale et al. (2008). 

different colour criteria, and with those of the SMGs, in the 
effort to elucidate the nature of the objects belonging to our 
candidate protocluster. 

3.2.1 Colour-colour diagrammes 

In the original paper by Magliocchetti et al. (2007) the dis- 
tribution of the Fsfiui/F24iiui versus the FR/F24^m colours 
has been shown as a powerful tool to discriminate between 
starburst dominated (^8Mm/^24(jm <1) from AGN domi- 
nated (-Fg Mm 7^24 Aim >1) objccts amoug the optical obscured 
24/im sources (see their Figure 2). The SED fitting results 
obtained in section 3.1 for all sources but J172030. 94-592310 
are consistent with the preliminary estimates obtained from 
the location of such sources in the [i^24-F(24/im)/i/RjR, 
i^24-F(24/im)/f8F(8^m)] colour-colour plane. Sources of 
this sample have \og{i'24F{24:^m) /{vRFn.)) > 1.74 and 
\og{v24,F{24:^m) /i^sF {8 fim)) > 0.55, corresponding to the 
colours of NGC 4418 (he avily obscured galaxy) shifted to 
a redshift of 2 (|Lutz et al. 2005), with the only exception 
of the AGNO which has log{iy24,F{2A^m)/iysF{8fj,m)) = 0.3 
and log{v24F{24:Hm) / (lynF^,)) > 1.74, putting it closer to 
the expected colours of the NGC1068 nucleus at a redshift 
around 2. 

The values of the 1.2mm/24/im flux ratios which are 
higher for the sources of the Spitzer sample of this work 
range from 3 to 5, except for the AGNO, which has 
F(1.2mm)/F(24/im) < 3. On a F(24^im) versus redshift 
plot these values coincide with those obtained by Lonsdale 



MAMBO Observations at 240GHz of optically obscured Spitzer sources 7 



Source name I fit Li,oi I fit II fit LsB II fit Lagn II fit 

(5 d.o.f.) Lq (4 d.o.f.) L0 Lg 



J172037.7+592359 


6.6 


6.5 X 10^2 


0.23 


4.1x10^2 


1.2x10^2 


J172030.9+592310 


81 


7.1 X 10^2 


23.9 


1.0x10^2 


O.SxlQi^ 


J172005.0+592430 


2.4 


5.7 X 10^2 


2.4 


3.4x10^2 


0.8x1012 


J171939.1+592512 


4.6 


5.3 X 10^2 


1.4 


2.5x10^2 


0.9x10^2 


J172025.2+592307 


2.5 


5.2 X 10^2 


1.1 


5.5x10^2 


0.8xl0i2 



Table 3. Results of the minimum x2 fits. Fit I is obtained adopting the Arp 220 SED (columns 2 and 3). Adding an obscured AGN 
component (fit II, columns 4-6) improves the values for three of the sources (we leave aside the source J172030. 9+592310 for which 
no acceptable fit was obtained; see text). We note that these results are only indicative as the number of photometric points is too small 
to allow any firm conclusion. The errorbars of the 24/im fluxes have been set to 20%. Errorbars on the total luminosities are estimated 
to be 30% - 35% (see text for detail). 



et al. (2008) for SWIRE sources (see their Fig. 8). Lonsdale 
et al. (2008) have observed a large part of the SWIRE sample 
at 1.2mm with the MAMBO array and have shown that the 
SWIRE sources and sub-mm selected galaxies (SMGs) have 
similar mm fluxes, but the latter have f (1.2mm)/f (24/im) 
ratios higher than the former by a factor of 2 to 10. Spitzer 
optically obscured sources have values intermediate between 
those of the SMGs and those characteristic of AGNs. This 
finding is independently confirmed by IRS spectra of SWIRE 
sources which show that part of the 24/im fiux is due to a 
hidden AGN. 

In Fig. H the location in the [f (24^m)/F(8/im), 
F(8/im)/F(4.5/im)] plane of our four sources with IRAC 
fiuxes are compared with those of the SMGs and of the dust 
obscured galaxies (DOGs; Pope et al. 2008). This colour- 
colour diagram was used by Ivison et al. (2004) as a di- 
agnostic tool to discriminate between starburst (SB) domi- 
nated and AGN dominated SMGs detected by Spitzer. Their 
SMGs turned out to lie in the locus of SB dominated galax- 
ies, except for two with AGN colours. The same diagram 
was used by Pope at al. (2008) to classify their DOGs 
{R- [24] > 14 and F(24Atm) > 100 ^tJy). These authors 
proposed the ratio F(8pim)/F(4.5/im) = 2 as the boundary 
between SB (lower values of the ratio) and AGN dominated 
SEDs. The conclusion was that 80% of the DOGs are star- 
burst dominated. 

All the sources analysed in this work but one (cfr Table 
1 and Figures O O have 8/im fiuxes below the sensitivity 
threshold of the relevant IRAC channel (i.e. 100 /iJy), yield- 
ing upper limits to the 8^m/4.5/im ratio in the range 2.2-3, 
and lower limits to the 24^m/8//m ratio in the range 3.8-4.4. 
For the measured 24/im fiuxes, this is indicative of a star- 
forming SED as an AGN/power-law SED most likely would 
have implied detectable 8/im fiuxes. In addition, adopting 
the Pope et al.'s (2008) criterion, our sources may have a 
SB component, though the 8/im observations are not deep 
enough to certain. The only exception is J172030. 9-1-592310 
which has 8/im/4.5^m ~ 3.6 and is therefore clearly in the 
realm of AGNs. 

IRAC photometry helps further in diagnosing the 
nature of these sources. In particular the [8/im/4.5/im, 
5.8/im/3.6/im] colours can be used to distinguish SEDs 
with evolved stellar population (which would determine 



a 'bump' in the SED between 3 and 8^m in the ob- 
served frame) and those with a power-law continuum all 
through the IRAC channels, characteristic of AGNs (i.e. 
Magiiocchetti et al., 2008a and Yun et al. 2008). Al- 
though we cannot draw any strong conclusion because of 
the lack of accurate IRAC photometry, we have at least 
upper limits on the flux ratios for the four sources de- 
tected by MAMBO, in the ranges logF(8/im)/F(4.5/.im) ^ 
0.3-0.5 and log F(;5.8/im)/F(3.6Mm) < 0.5-0.7, while 
the AGNO has log F(8^im)/F(4.5Mm) = 0.55 and 
log F{5.8fim) /F{S.6fim) ^ 0.2. The diagnostic diagrams of 
Yun et al. (2008) then confirm the conclusions above, i.e. the 
sources studied in this work are dominated by star forma- 
tion (SF) processes and have similar properties to the DOGs 
(Pope et al. 2008). 

As we have just seen, Spitzer data give strong hints on 
the nature of our sources. However, they do not allow us 
to unambiguously characterize them as they cannot provide 
direct information on their relationship with SMGs and on 
the importance of the AGN contribution to their SEDs. 
We can then go a step forward and exploit further infor- 
mation. For example, by using the mm and radio fluxes, 
together with the upper limits, we investigate the well- 
known, tight relationship between radio and far-IR emis- 
sions. The origin of this relation is attributed to star for- 
mation processes and has been extensively investigated 
for a quarter of a century. It is firmly established for 
the local Universe (i.e. Helou, Soifer, Rowan- Robinson 
1985) and appear s to also extend to higher redshifts 
(e.g. ICarrett 2002*; 'Ap pleton et al. 20041 : llbar et al. 20081 : 
ICarn fc Alexander 2009). 

By then comparing the [850/^m/1.4GHz, 850/im/24/im] 
colours of the Spitzer obscured sources with those SMGs 
having a radio and Spitzer counterpart we find that the 
latter have lower 850/.im/1.4GHz and 850/im/24/im ratios 
but both follow very tightly the correlation between the two 
colours, indicating once again that the sources in our sam- 
ple with measured MAMBO fiuxes are dominated by SF 
processes. 

The above correlation can also be exploited to estimate 
the redshift from the F(850^im)/F(1.4GHz) ratio, using 
the Carilh & Yun (1999) plot. We convert the MAMBO 
1.2 mm fluxes to 850/^m fluxes using the median ratio 



8 Paola Andreani et al. 



^850Aim/J'i.2mm — 2.5, as found by Greve et al. (2004), for 
SCUBA galaxies at z ^ 2. The source with 1.4 GHz radio 
flux has 850/im/1.4GHz ~ 22, placing it in the redshift 
interval 1-2.2. The 2 sources with radio upper limits have 
850/im/1.4GHz > 30 and > 45 and redshifts > 1.3 and 
> 1.5, respectively. 

As for the Spitzer source coinciding with the radio lobe, 
the convolution of the radio map with a Gaussian of FWHM 
corresponding to the 24/Ltm MIPS beam (6") gives a flux at 
1.4 GHz of ~ 3 mjy. The MIPS beam is comparable to the 
resolution of the radio VLA 1.4 GHz observations, which 
was 5". 

This value could be entirely associated to the lobe itself or 
may be due to the overlap of two components: one associated 
to the 24/im source and one to the lobe. We therefore con- 
sider this value as un upper limit to the real radio emission 
of this Spitzer source. In any case the correlation FIR-radio, 
which is the underlying assumption to infer the redshift, 
cannot be applied. 



4 DISCUSSION 

As is widely known, the flux of sources at a redshift around 2 
in the Spitzer 24/im channel may contain emission from the 
PAH feature at a rest frame wavelength of 7.7/im. Since this 
feature is associated to star forming processes, the sample 
selection based on the 24/im flux is biased as it enhances the 
probability of detecting star forming sources at this epoch. 

This flux band has been extensively used to estimate 
star formation in galaxies (see, e.g., Calzetti 2008 for a re- 
cent review and Rieke et al. 2009 for local galaxies). How- 
ever, the use of this band alone requires many a priori as- 
sumptions (e.g. comparison with low-redshift SEDs, PAH 
evolution, metallicity and the conditions in the photon dis- 
sociation regions) on the overall spectral energy distribution 
and completely misses star formation processes associated to 
molecular gas. We have estimated the star formation rate by 
fitting all the available photometric data on these sources. It 
follows that the uncertainties in this case are more related 
to the unavailability of the far-infrared fluxes, which can 
only be retrieved with future Herschel observations, and to 
the scant knowledge of the AGN contribution in the mid-IR 
domain. Because of their optical faintness {R >25.5) spec- 
troscopic redshifts of these objects are hard to obtain and 
they will be ideal candidates for ALMA follow-up observa- 
tions. 

Given the above difficulties, we tried to infer the na- 
ture of the sources in our sample from colour-colour plots, 
coming to the conclusion that the Spitzer optically obscured 
sources belonging to the analysed protocluster are likely to 
host an AGN and their mid-IR fiuxes may be affected by its 
emission. Quantifying this effect will be only possible when 
far-IR observations will allow us to derive the starburst con- 
tribution to the SEDs. 

One of the sources studied in this work has its position 
coinciding with one of the radio lobes of the radio loud AGN. 
At the locations of the radio core and of the second radio 
lobe there are no Spitzer detected sources. The convolution 
of the radio fiux coinciding with the location of the Spitzer 



source gives a fiux at 1.4GHz of ~ 3 mJy, which we consider 
as un upper limit to its real radio emission. 

The central source of this complex radio sys- 
tem is identified with FLSVLA172010. 2-^592425. 2 
jCondon et al. 2003h and at 8" from it there is an- 
other radio source, 7C1719-h5927 ijHales et al 1995l l. The 
fiux densities integrated over the whole system at 38, 330, 
610 MHz, and 1.4GHz are of 2000, 143, 33.7, and 35.6 
mJy, respectively. The 38, 330, and 1.4 GHz measurements 
are consistent with a power-law spectrum (5* oc i^"", 
with a ~ 1), while the 610 MHz flux density is a factor 
~ 2 below that. This atypical behaviour deserves further 
investigation. The central source has an optical counterpart 
with a magnitude of = 23.67 (Fadda et al. 2004). The 
optical redshift of this source is unknown. Using the R 
magnitude-redshift relationship for radio sources given by 
Rigby, Snellen, & Best [2007; their eqs. (9) and (10)] we flnd 
z ~ 1.55. Within the rather large uncertainties associated to 
this relationship, especially at z > 1 (see Fig. 6 of Rigby et 
al. 2007), this estimate is compatible with the radio source 
being at the estimated redshift of our Spitzer sources. If 
located a.t z = 1.73, the total luminosity at 1.4 GHz for 
the "concordance" cosmology (spatially flat universe with 
Qa ~ 0.7 and Ho — 70kms~^ Mpc~^), and a spectral index 
Q = 1, is 7 • 10^'^ ergs~^ Hz^^, i.e. a luminosity which is 
typical of an FR II (Fanaroff & Riley 1974) radio galaxy, 
consistent with its edge-brightened morphology. 

The total number of radio sources in the VLA-FLS area, 
as obtained by Condon et al (2003), brighter than 30 mJy is 
18. We then set, following Brookes et al. (2008), the fraction 
of radio sources at z>2 to be about 10%. By doing this we 
expect about 2 extended radio sources in the VLA-FLS field 
brighter than 30 mJy and with z>2. This number coincides 
with the two examples of extended radio sources associated 
to overdensities of Spitzer-obscured galaxies in the FLS (the 
other case is currently under study). We can then conclude 
that, despite of the small statistics there is a possible associ- 
ation of extended radio sources eX z ^ 2 with active mid-IR 
galaxies. 



5 CONCLUSIONS 

To test the possibility that the selected sources belong to a 
candidate protocluster and shed light on the nature of op- 
tically obscured galaxies defining the candidate filaments, 
we have carried out observations with the MAMBO array 
at the IRAM 30m antenna on Pico Veleta of five optically 
obscured sources. These latter trace one of the best exam- 
ples of a candidate filament, centered on a radio galaxy and 
roughly aligned with its jets. One of the Spitzer sources is 
located right at the edge of one of its radio lobes. Four of the 
observed sources were detected. Although the available data 
are insufficient to draw firm conclusions, the SEDs, obtained 
by combining our measurements with optical, infrared and 
radio observations, are consistent with the 4 detected sources 
being at the same z ~ 1.73 redshift and, therefore, trac- 
ing a real filament with a physical size of ~ 4Mpc, which 
may evolve into a rich galaxy cluster centered on a powerful 
radio-galaxy. 

The 4 detected sources most likely lie in the region of 
the [24/im/8/im, 8/im/4.5/im] colour-colour plot occupied 
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by dust obscured galaxies (DOGs; Pope ot al. 2008). The 
Arp220 SED provides reasonably acceptable fits for them. 
At the estimated redshift, the fits imply bolometric lumi- 
nosities in the range 5-7 • 10^^ L© and star formation rates 
of ~ 1000 Mq /yr. Their IR colours, however, arc bluer than 
those of sub-mm bright galaxies, suggesting either a signif- 
icant AGN contribution to IRAC fluxes or the presence of 
warmer dust, perhaps heated by young staxbursts. In fact, 
adding an obscured AGN component improves the x2 values 
for three of the 4 sources. The best fit AGN contributions 
to the bolometric luminosities range from 14 to 26% of the 
total. This is enough for the AGN to dominate the emission 
at 5.8, 8 and 24/im, while the SB dominates at 1.25mm and 
at A < 4.5/xm. The new fits arc again consistent with all 
these sources being at the same redshift, ~ 1.73. 

The source J172030.9+592310 may be a foreground ob- 
ject at z ~ 1.2, although the redshift estimate is very un- 
certain. Its SED is best interpreted as the combination of 
that of Arp220 with the "torus template" , corresponding to 
the model used to fit the SED of a heavily obscured type 2 
QSO, SWIRE J104409.95+585224.8 (Polletta et al. 2006). 
At z ~ 1.2, the best fit luminosities of the starburst and of 
the AGN components are 1.0 x 10^^ L© and 0.55 x 10^^ Lq, 
respectively. 

A large-scale interplay between AGN activity (such as 
that produced by radio jets) and star formation in the 
sources belonging to the protocluster may be at work here 
and may shape the properties of these sources. 
Following West (1994) we suggest that the protocluster may 
be in the process of forming the dominant cluster galaxy 
from the central powerful radio galaxy. 
To substantiate this rather speculative interpretation we are 
extending this investigation to other similar fields. 
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